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The transition from the microscopic Heisenberg model to the macroscopic elastic theory is carried 
out for the chiral magnetics of MnSi-type with the B2Q crystal structure. Both the exchange and 
Dzyaloshinskii-Moriya (DM) interactions are taking into account for the first, second, and third 
magnetic neighbors. The misorientations of the magnetic moments are divided into long-range and 
' short-range modes and it is shown that they are governed by different components of the DM vectors. 

The particular components of the DM vectors of bonds responsible (i) for the global magnetic twist 
and (ii) for canting between four different spin sublattices are found. It is demonstrated that the 
non-nearest magnetic shells could be very important because just the strong components of their 
£ ■ DM vectors contribute to the observed macroscopic twist contrary to the first neighbors. 
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C/3 ■ I. INTRODUCTION 



Chiral spin textures are studied now very actively for possible spin self-organization, unusual quantum transport 



phenomena and spintronic applications. A well established mechanism of spin chirality is the spin-orbit Dzyaloshinskii- 

rH ■ 

Moriya (DM) interaction which is responsible for intricate magnetic patterns in the MnSi-type crystals. Even half 
i ^ , ' a century after the discovery of the strange magnetic properties of MnSi [l|, |2j, the magnetics with the B20 crys- 
tal structure still amaze us with the variety and complexity of their magnetic phases and electronic properties 3j 
^ , contrasting the simplicity of their crystalline arrangement (only four magnetic atoms per a unit cell). Among the 
magnetic phases, both experimentally observed and hypothetical, are simple and cone helices [4|, [5fl , the Skyrmions 
: 1 ' ; ' 1 ■' . . . <• <>' - - 10- 12 similar to the blue 



and their lattices associated with the recently found A-phase [6H9J, possible 3D structures 



phases of liquid crystals, etc. This variety is due to, first, the lack of inverse and mirror symmetries, which gives rise 
^ ; to the chirality of the^stalline and spin structures; second, the frustrations resulting from the nontrivial topology 
of the trillium lattice [13| . that introduces a competition of various interactions between different pairs of atoms. 

Beginning from the discovery of the chiral magnetic properties of MnSi in 1976 0, Q till present day the most 
used approach to describe and predict twisted magnetic structures remains the phenomenological theory based on the 
Ginsburg-Landau free energy with an additional spin-orbit term first introduced by Dzyaloshinskii [llj, [ill, . 
However, the approach, which uses our knowledge of the system symmetry, is not able to say anything about the 
value of coefficients in the free energy, for instance, how they are connected with the real interactions in the crystal. 
The microscopic theories, e.g. the model of the classical Heisenberg ferromagnetics with a spin-orbit term originally 



developed by Moriya Uj, have in their turn the shortcoming that the number of variables, including the spins of all 



the magnetic atoms, is infinite, and therefore they are difficult to use for any analytical computations. Nevertheless, 
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even in spite of some doubts about the Heisenberg model in the itinerant magnetics, it is often used for digital 
simulations, if analytical methods based on the phenomenological theories seem to be too complicated [jjKjjj]. For 
the microscopic method it is of great importance to know the parameters of different interactions between pairs of 
spins: for exchange coupling of ith and jth atoms, and vectors Djj for DM interactions. Those parameters can 
be obtained in two ways: from comparison of theoretical results with an experimental data and by means of ab intio 
calculations. 

In our recent work fli] ]. using a coarse-grained approximation, the phenomenological constants J and T> of the 
simpliest elastic free energy have been connected with the corresponding parameters of the microscopic theory, which 
allows us to use the latter to obtain analitical results. A surprising detail was the recognition of inter-sublattice 
canting as a new microscopic feature of the magnetic structures in the MnSi-type crystals. Fig. Q] shows the difference 
between the twist and the canting by example of a ID spin chain. An experimental confirmation of the canting could 
give us the comprehensive answer on the question about the applicability of the Heisenberg model to the itinerant 
ferromagnetics of MnSi-type 2^]- But, even if the canting will not be found in the MnSi-type itinerant magnetics, 
there is a good chance to see it in dielectric chiral magnetics similar to recently studied multiferroics Cu2 0Se03 




22| , BiFe03 |23[ and langasite-type crystals [24 [ . 



Though the Heisenberg model with only nearest neighbors interactions is frequently used to simulate spin structures 
[12I fisL I19I ] , it does not always give an adequate physical description. Indeed, very often the coupling between second 
and next neighbors is comparable or even more considerable than that of the nearest neighbors. Thus, for example, 
in the weak ferromagnetic a-Fe2C>3, the antiferromagnetic exchange interactions and DM interactions are expected 
to be the most strong for third and fourth neighbors (see ab initio calculations in 25j); this is also supported by the 
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experimental data for the Heisenberg exchange parameters 

It also could happen that only with the interactions between non-nearest neighbors one can describe an experi- 
mentally observed phenomenon. For example, if in a ID spin chain with the ferromagnetic coupling between the 
nearest neighbors (Ji > 0) one turns on the antiferromagnetic interaction with the second neighbors (J2 < 0), 
then for IJ2I > Ji/4 a spiral magnetic structure appears owing to spontaneous breaking of the inverse symmetry. 
Such approach can explain the phase transition appearing at ^28 K and inducing ferroelectricity in the multiferroic 
Tb(Dy)Mn0 3 crystal 
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28|. The approach works also in the itinerant magnetics of MnSi-type. Just the competition 



between ferromagnetic coupling of the nearest spins ( J\ > 0) and antiferromagnetic coupling of the second and third 
neighbors (J2 ~ J3 < 0) was utilized to explain the observed alignment of magnetic helices along crystallographic 
directions (110) in MnSi at high pressure |29j. 

In this paper, we suggest a possible mechanism for effective reinforcement of the twist terms in the chiral spin 
structures of the B20 magnetics, taking into account the interactions with non-nearest neighbors. The reinforcement 
is induced by an interplay of the exchange coupling with the second and third neighbors and the canting [19I 
from DM interactions. In Sections [TT1 IVII the transition from the Heisenberg microscopic model to a continuous 
phenomenological one is performed, and the spin-orbit terms in the phenomenological energy are found up to the 
contributions of the order of (D/ J) 2 . In Sections I VIIHVIlTl the possibility of an extra twist induced by the canting is 
demonstrated. 
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II. EXCHANGE AND SPIN-ORBIT INTERACTIONS 

Both phenomenological and microscopical theories describing twisted MnSi magnetics contain terms induced by 
chirality of the system and associated with the Dzyaloshinskii-Moriya interaction of the spin-orbit nature. In the 
Heisenberg model, the extra term being associated with an individual bond, e.g. connecting magnetic atoms 1 and 2, 
can be expressed as 

Di2-[siXS 2 ], (1) 

where si and s 2 are the classical spins of the atoms, D12 is the DM vector of the bond 1-2. In principle, there could 
be magnetic moments at Si atoms but this effect will be neglected below. Let us shortly describe the properties 
of the DM vector [17j . (i) As it is obvious from Eq. ([1]), the sign of the DM vector depends on which atom of the 
bond we consider as the first. Indeed, because the cross product changes its sign when Si and s 2 are rearranged, then 
D21 = — D12. (ii) The structure changes its chirality under inversion, whereas the energy remains unchanged, which 
means that D12 is a pseudovector. (hi) The DM vector of a bond possesses the local symmetry of the bond. Thus, 
in MnSi the DM vectors are of their most general form, i.e. they have three independent components, because the 
Mn-Mn bonds in the B20 structure do not possess any internal symmetry element, (iv) The DM vectors vary from 
bond to bond; the DM vectors of the equivalent bonds have the same length but may be of different orientation (just 
the case of MnSi). In particular, if two bonds are connected by the rotation symmetry transformation of the space 
group P2i3 of the crystal, then their DM vectors do by the corresponding rotation of the point group 23. 
In the phenomenological theory, the chiral interaction is induced by the following extra term 

VM • [V x M] (2) 

in the expression for the magnetic energy density. Here M is a continuous field of the magnetic moment, T> is a 
pseudoscalar constant of the interaction. 

Because both the theories describe the same matter, there should be a relationship between them. In particular, 
the terms ([I]) and (J2|) of the different approaches should be somehow connected. Indeed, it was shown in Ref. [ijj 
that in the nearest neighbors approximation the constant T> of the phenomenological theory is proportional to the 
component (D x — 2D y — D z ) of the DM vector of the bond (—2x, \-,\ — 2x) between neighboring m ang anese atoms. 
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34], the DM 



However, there is a problem here. According to different spin-orbit calculation schemes 
vector associated with a bond should be perpendicular or almost perpendicular to the bond. In the present case 
it means that the component {D x — 2D y — D z ), which lies almost along the bond, constitutes only a small part of 
the DM vector length. Taking into account that having the relativistic nature spin-orbit DM interaction serves as 
a small additive to the ferromagnetic exchange coupling, we can conclude that the twist observed in the MnSi-type 
crystals, particularly in MnGe, seems to be abnormally strong. A possible solution of this problem is that the nearest 
neighbors approximation is not sufficient and hereinafter in the paper we develop this idea in details. 

In the following sections we will show, how to perform transformation from the microscopic Heisenberg model to a 
macroscopic elastic theory. 



4 



III. FROM THE HEISENBERG MODEL TO A CONTINUOUS APPROXIMATION 

In the classical Hcisenberg model of magnetics with an extra interaction of the Dzyaloshinskii-Moriya type, the 
energy of the system is expressed as a sum of pair interactions between magnetic atoms and the interaction of 
individual atoms with an external magnetic field H: 

E = E ( E {" J ^ ■ s i + D ».« • k x s ^ - ^'tr^'^ l - ^bH ■ £ * ] . (3) 

cells \n,{ij} n ) j=1 J 

Here the external summation is over all the unit cells of the crystal, n enumerates magnetic coordination shells, 
the sum {ij} is over all the bonds of the shell n, the sum over i is taken over all magnetic atoms in the unit cell; 
Jn = Jn — (D^/AJn), J n is the exchange coupling interaction parameter of the nth shell, T)n,ij is the DM vector of 
the bond ij of the nth shell, Sj are classical spins of magnetic atoms, |s| = 1, g[is is an effective magnetic moment of 
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the atom. The terms of the order of D\ are usually ignored but sometimes they can be important 

The MnSi-type crystal has £?20 structure (P2i3 space group) with four magnetic (Mn) atoms occupying crystallo- 
graphically equivalent 4a positions within a unit cell, ri = (x,x,x), r 2 = (1 — x, \ + x, ^ — x), r 3 = (A — x, 1 — x, h + x), 
r4 = (i + x, h — x, 1 — x) 37]. The shortest bonds between the magnetic atoms have the length J h — 2x + 8x 2 . 



The next environment consists of two close magnetic shells with the radii ^/ § — 6x + 8x 2 and \J j + 2x + 8x 2 . The 
parameter x defines which of them is closer to the initial atom. When x > i (a; « 0.138 in the case of MnSi), then 



the shell with the radius \/| — 6x + 8x 2 is closer, so we will refer to it as to a second one. When x = r g J the 2nd and 

38j, connecting 



3rd neighbors are at the same distance so that the manganese sublattice gains the space group P4332 
the 2nd and 3rd shells by a symmetry transformation. This could establish a linkage between the DM vectors of the 
shells, but the silicon sublattice (having in its turn the space group P4i32, when irs^ = |) breaks the symmetry and 
the DM vectors of the shells are different even in this case. If x < |, then the 3rd neighbors become closer than the 
2nd ones. 

Every magnetic atom possesses six neighbors in each of its three nearest magnetic shells, so there are 36 bonds in a 
unit cell (in twelves for each kind). Fig. [2] shows the magnetic environment of a manganese atom in MnSi. The atom 
lies on a 3-fold axis. Each of the first three magnetic shells consists of two equilateral triangles. All the atoms of the 
2nd shell are on one side of the plane perpendicular to the 3-fold axis, all the atoms of the 3rd shell are situated on 
the other side of the plane. 

Let (Di x , Di y , Di z ) be the coordinates of the vector D112 corresponding to the bond bi.12 = (— 2a;, i, i — 2x) 
directed from ri to r 2 — (1,0,0), (£>2x, D2y, -D22) be the coordinates of the vector D242 corresponding to the bond 
b2,i2 = (1 — 2a;, i, i — 2a) directed from ri to r 2 , and (D 3x , D 3y , D 3z ) be the coordinates of the vector D 3i i 2 
corresponding to the bond b3 i i2 = (— 2x, ^, — \ — 2a) directed from to r 2 — (1, 0, 1). 

We can specify directions for all twelve bonds of a shell in the way that all directed distances by between neighboring 
atoms will connect to each other by the symmetry transformations of the point group 23, so corresponding DM vectors 
will be connected by the same transformations. In the MnSi-type crystals there is only one type of magnetic bonds 
directed from a central atom of the type t to an atom of the type t' 7^ t at a given shell (n = 1, 2, 3). We can take 
advantage from this fact and introduce for each pair ij (i, j = 1, 2, 3, 4, i ^ j) a local triad 



(n - Tj) _l (n + Tj) ± [n x -n,-], 



(4) 
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and then use it as a basis for vectors by and Dy . Here Ti is a vector directed along the 3- fold axis passing through 
the position iv 

Tl = (1,1,1), 

T 2 = (-1,1,-1), 

V ; ' (5) 

T 3 = (-1,-1,1), 
T 4 = (1,-1,-1). 

Then the bond vector and the DM vector of an arbitrary bond in three first magnetic shells can be written as 

bi,y = (1/8 - x)(Ti - Tj) + {n + Tj)/A + [t< x Tj]/8, (6a) 
b 2 ,y = (3/8 - x)(Ti - Tj) + {n + Tj)/A - x Tj]/8, (6b) 
b 3 ,ij = (-1/8 - x)(75 - 7j) + + r,)/4 - [75 x Tj -]/8, (6c) 

D, w - (r* - t 3 ) + ^(T i + Ti ) -^[tjx r,], (7) 

with Z3 n ± = D na; ± -D„ 2 ; here D„^j 7^ ~D n ^j because the indices ij and ji designate two different bonds. 

Notice that the energy ([3]) does not depend on the parameter x, so it can be chosen arbitrarily, making convenient 
transition to the continuous approximation (Fig. [3Jb)). Indeed, to provide the transition the continuous unimodular 
vector functions will be defined, §1, §2, §3, S4, which coincide in special points with the spin values in corresponding 
atomic positions. If the parameter x is equal to its experimental value, e.g. x — 0.138 in MnSi, then we should take 
the functions in the real atomic positions in order to obtain the real spins. In the case of an arbitrary choice of x we 
should take the function values in the points shifted from the real atomic positions. 

For our calculations it is convenient to choose the parameter x from the condition 

Mn -Tj)® b ntij = 0. (8) 

n,{ij} 

Notice that this 3 x 3-tensor is proportional to the unit one due to the averaging over the symmetry elements of the 
point group 23 (see Appendix, Eq. (|A4[)). 

Thus, for example, in the nearest neighbors approximation the condition gives x = |, which is shown in Ref. [3| 
to be necessary in order to obtain smooth spin functions. 

For n = 1, 2, 3 using Eqs. (|A"8|) we find from Eq. © 

1 w 

*=8 + 4' (9) 

W = — = — . (10) 

J1+J2+J3 

Notice that the parameter w does not depend on the freedom in the choice of x, because it is expressed only through 
the interaction constants. Hereinafter we suppose that \w\ < 1. 

The bond ij of the nth shell connects the function values Sj(r) and Sj(r + b Hi ij). Using the continuity of the 
functions s we can write 



s,-(r + b„,y) = exp(b„,jj • V)sj(r) 



(11) 



6 



and pass from the summation over the unit cells to the integration over the crystal volume (in chosen units the lattice 
parameter a = 1 and the unit cell volume is equal to 1): 

E= J v dr (j2 n ,{ij}{-JnSi ■ exp(b n ,y • V)sj + H n>i j ■ [s 4 x exp(b IW ■ V)§j] 

4 \ (12 J 

-(Pn,ij ■ &i)(Pn,ij ■ exp(b niij • V)Sj)/(2J n )} - j/i B H • £ i=1 M • 

IV. MAGNETIC MOMENT DENSITY, CALCULUS OF VARIATIONS AND PERTURBATION 

THEORY 

In the phenomenological theory, the magnetic moment density M is used as an order parameter and considered as 
an experimentally observed physical quantity. Nevertheless, when we try to define a smooth continuous function M 
from a discrete distribution of spins of the magnetic atoms, an ambiguity arises from the fact that we can determine 
the weight function in different ways when averaging the spins on a local volume. In this case, the ambiguity results 
in that the smooth functions S; having specified values in a discrete set of points can be defined by infinite number of 
ways (Fig. [Ha)). We will avoid the problem, supposing that the functions Sj are defined in the most convenient way. 

It would be natural to define the average magnetic moment per atom as 

M = 5/i B m, m = + s 2 + s 3 + s 4 ). (13) 
Let us introduce the vector functions u x , u v , u z defined by the relations 

U x = (Si - S 2 - S 3 + S 4 )/4, 

Uj, = (si + s 2 - s 3 - s 4 )/4, (14) 
u z = (§i - s 2 + s 3 - s 4 )/4, 

so that 

ha — m a + T ia u aa , a,a = x,y,z. (15) 

Hereinafter the summation on repeated indices is implied. Four conditions |sj| — l,i — 1,2,3,4, can be rewritten 
using the invariants 

io = 1(4 + 4 + H + si) 

Ix = — s\ — s\ + sf) 
ly = + *2 — s l — *I) 

Iz = — *2 + *3 — *1) 

Let us rewrite Eq. (fT2|) using just defined functions and include in the integrand the Lagrange terms Aoio, X x Ix, 

Ay-^y, X Z I Z . 

E = f v dr{ S„,{y}h^»( TO « + T icrU aa ) exp(b n ,y ■ V)(m Q + T jp U pa ) 

(mp + n a u a p) exp(b n ,ij • V)(m 7 + t jp u p1 ) 

)D n ,ijp exp(b n ,y • V)(mp + T jp u pf} )] (17) 
-4#//bH • m + A (m 2 + u aa u aa ) + X x (m ■ u x + u y ■ u z ) 
+X y (m ■ u y + u z ■ u x ) + X z (m ■ u z + u x ■ u y )}. 



m + u x 
m • u r - 



m • u 



m u. 



0. 
0. 



UX ■ Uy = 0. 



(16) 
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In order to express the energy as a functional of the magnetic moment m, E should be minimized by u x , u y and 
u z using the calculus of variations. The variation of the energy (| 1 7[) is 



SE = / 5£dr = / dn aa ^ aa dr, (18) 
Jv Jv 

^<ja = J2n,{ij}{-Jn[Tio exp(b„,y • V) + Tj a exp(-b n> ,j ■ V)]TO Q 

-J n [n<7Tj p exp(b bn.ij • V)]zi pa 

Sa^DnjjplTia exp(b n ,y • V) - Tj CT exp(-b n ,,j • V)]m 7 
-ea/37-£ > n,ij/9[T«r7j>exp(b rl) ij ■ V) - T icr T lp exp(-b„,ij • V)]u P7 

Dn^jp lTia exp(b n ,ij • V) + Tjo- exp(— b^ij • V)]m^ (19) 

V )K/?} 

^zU xa + \ x ii za , cr = y, 

Here the evident equation was used, which corresponds to an integration in infinite volume, 

J /(r) exp(b • V)g(r)dr = J g(r) exp(-b ■ V)/(r)dr. (20) 

The minimum of the energy is determined from the condition (5£ = with arbitrary functions 8u aa , so the problem 
is reduced to the system of 9 equations f m — 0, a, a — x,y,z, with the extra conditions (|16p determining the 
functions A (r), A^r), A y (r), A z (r). 

The general solution of the problem is too difficult, if possible, but we can use the perturbation theory with a small 
parameter 

D/J<€1, (21) 

which is the ratio of typical values of the spin-orbit D and exchange J interactions. Let us assume that this parameter 
also describes the typical order of magnitude of derivatives on coordinates and the spin components responsible for a 
canting: 

|V|~|iv|~D/J. (22) 
Let us use the consecutive approximations in order to find a solution in the form 

u CT = u«+u( 2 >+u( 3 >..., 

X = X^+X^+xfK.., (23) 

A CT = Act ' + Act ' + Act • • • , 

where upper index (p) means that the corresponding term is of the order of (D/ J) p . 



V. THE FIRST APPROXIMATION 



The zeroth order equations on D/ J have the view 



M T ™ + T ja )m a = (24) 

n,{ij} 
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and become trivial after the summation on bonds, see Eq. (|A6a ). Let us write down the system ^ aa = keeping 
only the 1st order terms on D/J: 

/^n,{ij} i Jn(j~ia ^~ja^n,ij ll^ct Jnij'ia^'jp Tj&T-ip^Upa (25) 
-Ea^D^ij^Tia ~ T jcr )m 7 } + 2X^Ual + \a m a = 0. 

The first summand in curly brackets is zero in accordance with Eq. ([5]), two other can be calculated using Eqs. (1A7F ) 
and CM)- Then, 

2(A[, 0) + 4( jj + J 2 + J 3 ))u« + 8(D 1+ + D 2+ + D 3+ )e aaimi + \V>m a = 0. (26) 
The normalization conditions are 

"M = l-™ 2 , (27) 



= 0, (28) 



where the first equation is of the second order on D/J. 

Multiplication of Eq. ([26]) by m a and summation on a with use of Eq. ([28]) give 



A« = 0, (29) 

therefore 

m _ 4(D 1+ + D 2+ + D 3+ )e 

A[, 0) +4(J 1 + J 2 + J 3 ) 

The substitution of Eq. d30j) into Eq. ([27]) gives 

^ + 4(Ji + J 2 + J,) = ^i + + ^ + + ^| m> (31) 



\fl — m* 



and, finally, 
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u^a = ~sign(D 1+ + D 2 + + D 3+ ) 1 m £ga7 m 7 , (32) 

V2m 



where the sign of the right parts of Eqs. pip and (|32[) is chosen to be seen in future from the condition of the minimum 
of the canting energy. 

The substitution of Eq. ([32]) into Eq. ([15]) gives 



Sj = m + n[Ti x m], (33) 

with 



K = sign(D 1+ +D 2+ +D 3+ )^-= . (34) 

V2m 
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VI. ENERGY DENSITY 

Let us begin from the calculation of the contribution to the energy density from the magnetic moments m and their 
derivatives: 

£<W = - J n m a m a = -12(Ji + J 2 + J 3 )m 2 , (35) 
n,{ij} 



£m = X! {~^« TO a( b 'w ' + e a p 1 D n . i j OL mf i 'm~ l ) = 0, (36) 



£m } = X) n ,{ij}{-5^n m a( b «,v) ' V) 2 m a + £ Q/ 3 7 Al,ija TO/3 (b„,ij ■ V)m^ 

D n ^ 3 pm a mp}. 
Using Eqs. (|A11|) . (|A13|) . and (|A10|) we obtain 

4? } = JmAm + £>m[V x m] - — ± + — ^ + — ^ m 2 , 

\ Ji J 2 J3 / 



(37) 



(38) 



with 



j _ 3J 2 + 3J| + 3Jf + lOJi J 2 + lOJi J 3 + 22J 2 J 3 
4(Ji + J 2 + J 3 ) 



(40) 



V = 2w(D 1+ + D 2+ + D 3+ ) - 4Di,ia ■ bi iia - 4D 2 , 12 • b 2il2 - 4D 3 , 12 • b^ 12 

= 2w(D 1+ + D 2+ + D 3+ ) - (-£>!_ + 2Di v ) - (2D 2+ + D 2 _ + 2D 2y ) - (-2D 3+ + D 3 _ + 2D 3y ). 

Notice, that the first term changes its sign simultaneously with J 2 — J 3 . 
The first term in Eq. (1381) can be rewritten as 

orp orp 2 

where the second term gives a contribution to the surface energy only: 

-\j I drV ■{Vm 2 ) = -\j f, df • Vm 2 . (42) 
2 Jv 2 J s 

Far from the transition between the paramagnetic and ferromagnetic states, the absolute value m of the magnetic 

moment changes slowly, and this contribution to the energy can be neglected. However, near the phase transition, m 

can undergo considerable changes, and in the crystals with a significant surface (nano crystals, thin films) the term 

(|42p could play an important role. In par ticular, it could be important for the stabilization of the A-phase observed 

in thin films of Feo.5Coo.5Si and FeGe [3^, Q- 

The first term in Eq. (jH]) defines the deformation energy of ferromagnetic ordering, when J > 0. When J < 0, the 

continuous approximation is not applicable any more owing to the antiferromagnetic tendency of spins. The criterion 

J > of the applicability of the continuous approximation can be represented geometrically (Fig. |4|) using the formula 



7 Zj ± + %k±y±Z&: (43) 

J -l Jl TTT+ ' ( } 
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with j± = (J 2 ± J 3 )/Ji. 

Let us find now the contributions from the canting into the energy density: 



£ u ] = - E J «( r - + r *>£l m « = 0, (44) 



= 4(Jj + J 2 + J 3 )uScUla + 8(L>i+ + D 2+ + £> 3+ )e (T; 3 7 M^ /3 ) m 7 
= 4(Ji + J 2 + J 3 )(l - to 2 ) - 8V2|£>i+ + L» 2+ + %|Wl-m 2 . 



(45) 



Thus, we can finally rewrite the bulk energy density as a function of the magnetic moment m accurate within the 
second order terms on D/J: 



£ = ^TFTT^f + x m] + {.h + J 2 + J 3 )(4 - 16to 2 ) 

'2D? , 2D? 2rt\ o . „ ( 4b J 



iV^|Di+ + D 2+ + D 3+ \mVl-m 2 - (H^i + ^i + ^£) TO 2„ 4 5MsH • m. 



Let us consider the expression in details. First two terms with derivatives are nothing but the deformation energy 
of the conventional phenomenological theory of the chiral magnetics. However here the values J and T> are expressed 
through the parameters of the microscopic theory accordingly to Eqs. (|39j) and (|40[) . The following term, on conditions 
that to ~ 1, gives — 12 (Ji + J 2 + J3), that is the energy of 36 ferromagnetic bonds in first three shells in the unit 



cell. Then, the term with vl — to 2 is a contribution of the canting. It is always negative, which means that the 
canting is an important and unavoidable peculiarity of the magnetic structures of the MnSi-type crystals. When 



minimizing the energy on to, this term gives the contribution to the derivative proportional to 1/Vl — m 2 , which 
becomes dominating, when to — > 1 — 0, impeding to to exceed 1. 



VII. EXTRA TWIST INDUCED BY CANTING 

Let us return to Eq. (jTOj) defining the DM parameter T> of the phenomenological theory. As it was mentioned above, 
from the physical point of view, the DM vector is expected to be almost perpendicular to the corresponding bond 
and, consequently, last three terms are expected to be small. Some speculative estimation can be made with the use 

based on the Moriya theory [l7|, which gives us following expression for the DM 
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of the well known Keffer rule 
vector 

Di3=%xr 2i ], (47) 

where D is an unknown coefficient, and the vectors ru and v 2 i are directed from the positions of 1st and 2nd magnetic 
(Mn) atoms to that of the intermediate nonmagnetic (Si) atom realizing the spin-orbit interaction. It is evident that 
D 12 -L r i2 = r 2i - r u . 

If the structure parameter x = i, then b„ = h' n and the components — D\- + 2Di y , 2D 2+ + D 2 - + 2D 2y , and 
— 2Z?3 + + D3- + 2Dz y are exactly parallel to the corresponding bonds. In the known crystals the parameter x is close 
to |, e.g. in MnSi x — 0.138 and the angle between the direction [121] corresponding to the component — D\- + 2Di y 
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and the direction of the bond b>i 12 = (—2a;, 2a;) is about 3.4°. The more particular conclusion from the Keffer 

rule is that [2^] 

D Ma ~(2-4ic,l,4ar-l) (48) 

and, thereby, the projection of the DM vector on the direction [121] makes up only 2.4% of its length, so we can 
conclude that 

- D x _ + 2D ly ps 0. (49) 

The similar estimations can be made for the DM vectors of the 2nd and 3rd magnetic shells. Thus, the manganese 
atoms in the positions i"i and T2, connected by the bond b2,i2) have two common silicon neighbors in the positions 
(| — Xsi, 1 — x$i,xsi — 5) and {xgi — 5, § —x$i, 1 — xsi) {x$i — 0.845) and form together with them a slightly distorted 
parallelogram resulting in T>2,i2 ~ 0. On the other hand, the manganese atoms in the positions ri and T2 — (1, 0, 1), 
connected by the bond h 3 ^2, have only one common silicon neighbor (1 — xsi, xsi — hi \ ~ xsi), which results in 

D 3 , 12 ~(-4s 1 -l,4a:-l) (50) 

(in the units of Eq. (j48|l ). so ID3 1 / 1 Di | « 0.67, and the projection of the DM vector 03,12 on the direction [123] makes 
up only 2.2% of its length. 

Taking this into account and neglecting the products D n • , the constant T> responsible for the twist of the chiral 
magnetic structure can be approximately written as 

v „ 2(J 2 ~ J 3 )(D 1+ + D 2+ + D 3+ ) (51) 
J1 + J2 + J3 

The combination D\ + + Z3 2 + + D 3+ in the numerator of the expression is responsible for the canting in accordance 
with Eq. pop , and J2 — J 3 is a combination of the exchange interaction constants between non- nearest neighbors. 
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2, 



Therefore, the canting, initially considered as a supplementary microscopic peculiarity of the chiral magnetics 
can be in fact an essential cause of the extra twist experimentally observed; in the next section we will consider a 
speculative example clarifying the situation. 

This fact can be reflected in a relationship between observable values. Indeed, in the absence of the magnetic held 
the minimum of the energy corresponds to the helix with the propagation number 

k = — ^ 4(J 2 - J 3 )(A+ + £>2+ + £>3+) m 

2 J ~ 3 J\ + 3 J 2 2 + 3 J 3 2 + lOJi J 2 + 10 Ji J 3 + 22 J 2 J 3 ' 

k = ±|k|, where k is the helix propagation vector, the sign defines the chirality ("plus" for the right-handed helix). 

In the case of the strong magnetic field (H > i?2c) the helix becomes unwound, so m = const, m || H. The 

minimization of Eq. (|46j) on m gives 

VT^ « ^Di+ + D> + + D, + \ (53) 
i(J 1 + J 2 + J 3 +9ti B H/8) 

The substitution of Eq. ((53]) into Eq. ((34)) gives 

D 1+ + D 2+ + D 3+ 3 J 2 + 3J| + 3J3 2 + lOJi ~h + 10Ji J 3 + 22J 2 J; 



3 k. (54) 



4( Ji + J2 + J3+ g^ B H/8) 16( J 2 - J 3 ) ( Ji + -h + J3 + gu B H/8) 

Thus, related are the propagation number k of the helix magnetic structure and the canting constant k in the 
unwound state. The latter can be mesuared, for example, in neutron or X-ray diffraction experiments by studying of 
the forbidden reflections 00£(£ = 2n + 1) [2oj|. 
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VIII. SPECULATIVE DISCUSSION 

Let us demostrate by a simple albeit not very realistic example how the interactions with non-nearest coordination 
shells can result in an essential twist gain. Let us consider a periodical ID chain of spins with a local interaction 
between them, composed of unit cells containing two spins, say 1 and 2 (Fig. [Ha)). Let us assume that all the spins 
can rotate in a plane, and write the energy of the chain as a function of differences of the angles, 



E = J2 



{ Cl {^ ] - ^ - P) 2 + Orf - A - Pf + M<pP - ^l) 2 + MvP ~ &) 2 } ■ (55) 



Here cp^ and ipffi are the orientation angles of the spins 1 and 2 of the nth cell, laid off from an arbirary direction; 
C\ and Ci are positive constants. The condition C\ ^ Ci determines the chirality of the structure. 

When J2 = 0, Eq. ([55"T) describes a periodical magnetic structure with the angle (3 of the canting between the 
spin sublattices 1 and 2 (Fig. [5{b)). We suppose that this structure is a result of the competition of two ordering 
interactions between neighboring spins, the ferromagnetic one and the twisting chiral one; besides, there is a reason, 
which is not considered here, eliminating the twist in this state. 

The question arises: what happens with the structure after introducing of ferromagnetic (J2 > 0) or antiferro- 
magnetic ( J2 < 0) interaction with the second neighbors from the complementary sublattice? The minimization of 
Eq. (|55|) gives the solution in the helix form 



(!) X 

fn = no, 

( 2 ) X 

iph. — no - a, 



(56) 



where 



c l + c 2 C2 + J2 § (57) 



C1 + C2 + 2J2' C1 + C2 + 2J2 
is the tilt angle between the spins in the unit cell, and 

x Hc 2 -Ci)p 

9J 2 2 + 5J 2 (C 1 +C 2 ) + CiC2 V ' 

is the twist angle per one period of the chain (Fig. (He)). 

At the first sight it seems to be paradoxial that the introducing of an additional aligning interaction induces a 
twist, but a simple analysis of the fact shows that there is no contradiction here. When IJ2I — ► 00, the angles a and 
5 go to zero as expected. In order to understand the system behavior for the finite values of J2, let us consider the 
numerator of Eq. ([55)) . The factor (C2 — C\) reflects the degree of the internal chirality of the structure. The product 
J2/3 reminds of the lever torque, with /3 playing the role of the lever arm and J2 being analogous to the rotating force. 
We can imagine that the "aligning force" J2, being applied to the initially tilted by the angle j3 sublattices, results in 
structure distortions, which in their turn induce the twist due to the potential chirality of the structure (C\ ^ C2). 
Note, that the change of the sign of the constant J2, corresponding to the transition from ferromagnetic (J2 > 0) 
to antiferromagnetic ( J 2 < 0) coupling of the second neighbors, increases the degree of the twist (reduces the helix 
pitch) and changes its handedness. 

We expect that the similar effect takes place in the chiral magnetic MnSi. Indeed, there is an evident similarity of 
Eqs. (fS"Tj) and {55J. In this case the role of a lever is played by the canting induced by the DM interactions, whereas 
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the "force" is the anti- or ferromagnetic aligning interaction with the atoms of non-nearest magnetic coordination 
shells. It should be noted that three first magnetic shells in the B20 structure contain the atoms belonging to the 
sublatticcs different of that of the central atom (Fig. [5]). 



IX. SUMMARY AND DISCUSSION 

An essential difference of the microscopic description of the magnetic properties of the MnSi-type crystals from 
the phenomenological one is the usage of the pseudovector D instead of pseudoscalar T> when describing the 
Dzyaloshinskii-Moriya interaction. The presence of the extra parameters (pseudovector components) results in the 
existence of the local canting, the feature not studied yet in the magnetic twisted structures. An important problem 
solved in the present work is how to distinguish the components of D responsible for the twist and the canting. In 
Ref. [18[ it was shown numerically (in the nearest neighbors approximation) that responsible for the twist were the 
components of the D-vectors lying along the bonds, whereas for the canting did the D components directed perpen- 
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dicular to the bonds and lying in the planes of bond triangles in the trillium lattice. In Ref. 
and showed that the real components of D inducing the twist and the canting lay along crystallographic directions 
closed to those found in Ref. 18[- However, the fact that the D components responsible for the twist lay almost along 
the bonds could still seem accidental. Here we show that it is also true for the 2nd and 3rd neighbors. Indeed, it 
follows from Eq. (|Al2j) for the small parameter w that the twist is defined by the scalar product of the DM vector D 
and the bond b', with the additional condition x = g (in real MnSi-type crystals x « g). 

Generally speaking, for symmetry reasons the DM vectors of first magnetic shells in MnSi have to be of general 
view, and we can say nothing about how the extents of the twist and the canting can be connected to each other. 
Nevertheless, it is expected from the quantum mechanical description of the Dzyaloshinskii-Moriya interaction that 
the DM vectors are almost perpendicular to the bonds. This means that in point of fact the products D • b' are tiny, 
and, even in the case of a small w, the twist can be defined by the components D + of DM vectors (see Eq. (|A12|) ). 
which are due to Eq. (|A9p to be also responsible for the canting. It leads to the conclusion that the canting, initially 
being considered as an additional microscopic effect in relation to the global twist, can in fact serve as a cause of the 
abnormal twist experimentally found in the MnSi-type crystals, particularly in MnGe. It is also important that the 
contribution of non-nearest magnetic neighbors should be taken into account. 

In the simple phenomenological theory describing twisted magnetic structures, it is supposed that |M| = const, 
which is roughly true at low temperatures, far below the phase transition from the paramagnetic state to the twisted 
ferromagnetic one. However, this condition makes energetically unfavorable such structures as the Skyrmions and 
their lattices, associated with the A-phase observed close by the transition point. In order to overcome this problem as 
well as to describe the critical phenomena, two additional terms, M 2 and M 4 , limiting the value of magnetic moment 
are included into the free energy [l^. The presence of the terms decreases M in the regions with a large density of 



the magnetic energy, thereby decreasing the energy of the whole structure. Thus, in Ref. [42j the terms M 2 and M 4 
are used in order to calculate the energy of Skyrmions. It is found that M decreases considerably nearby the core 
of the Skyrmion, where the magnetic moment M has the opposite direction to the external magnetic field, and has 
a maximum at some distance from the core, where the energy gain from the double twist is maximal. However, in 
the latter case the magnetic moment exceeds its saturation value Mo, which is not acceptable for physical reasons. 
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In the present work we show that there should be a contribution from the canting ~ m\J\ — m 2 into the energy 
density, which can play the same role as M 2 and M , but does not allow the magnetic moment modulus to exceed 
the saturation threshold. Besides the canting, the thermal oscillations of spins also contribute to the reduction of 
the magnetic moment M. Nearby the transition point the amplitude of the oscillations can be comparable with the 
canting. Moreover, the less is the effective local field it, e ff,i — —dE/dsi, acting on the individual spin s,-, the more 
oscillations and cantings. Therefore, the reinforcement both of the thermal oscillations and the canting have the same 
cause, so they should give a similar effect. 

If the canting would be observed in MnSi, the direct confirmation of non-nearest neighbors effect could be possible 
using Eq. (|54[) . connecting the propagation number of the magnetic helices with the magnitude of the residual canting 
in the unwound state in a strong magnetic field. The equation involves only the exchange constants J and does 
not depend on the DM vectors. The possibility of an experimental proof of the theory should stimulate ab intio 
calculations of the interaction constants. In particular, the ab intio calculations could give the precise estimation of 
the difference J 2 — J3 defining the size of the parameter w responsible for the extent of the twist. For example, the 



authors of Ref . 



29| suppose that J2 « J3 due to the closeness of the 2nd and 3rd magnetic shells (in MnSi x — 0.138 



and 63/62 ~ 1.06). If | J2 — J3I <C 1, then the contribution of the 2nd and 3rd neighbors can be insufficient to explain 
the strong twist observed in the MnSi-type crystals, especially in MnGe. Here the realistic calculations would give a 
tip about the direction of the further search. 

It has been found in Ref. 43] that the propagation number k of the magnetic helix in the alloy Fei-^Co^Si is 
strongly dependent on the concentration x of the cobalt atoms. Thus, when x changes from 0.05 to 0.15, the helix 
period decreases abruptly in several times. In the present work we show that k is proportional to the difference J2 — J3 
of the exchange interactions with the atoms of the 2nd and 3rd magnetic shells, and J2 ~ J3 [29]. We suppose that 
the considerable increase of spiraling could be induced, for example, by a different filling of the cobalt atoms in the 
shells. This hypothesis can be tested experimentally. 

Another interesting fact, not being yet explained within the framework of microscopic theories, is the helix ordering 
along some special directions, e.g. (Ill) or (100). Usually this ordering is associated with a weak anisotropic exchange 
[15) 1 . but in fact the ordinary isotropic DM interaction also can result in the appearing of anisotropic terms in the 
energy of spiral orientation in the lattice with a cubic space group. These contributions of the order (D/J) 4 will 
define the critical magnetic field H&, at which the helix comes off from its preferable zero- field direction. This is 
in a good agreement with the observed ratio of the first and second critical fields H c i <C H C 2, because it follows 
from the theory that H c i/H C 2 ~ (D/J) 2 . Indeed, e.g. the period of the magnetic helix in MnSi makes about 40 
unit cell parameters, which gives the value 27r/40 for the propagation number modulus |fc|. On the other part, it 
follows from the phenomenological description that k = T)/(2J), or D/J ~ 7r/10 ~ 0.3. This gives us the estimation 
Hci/ 1 H C 2 ~ 0.1, which is in a good agreement with the experimental data. In our previous works [lill^ we proposed a 
coarse grain approximation, which allowed to calculate only the contributions to the energy proportional to (D/J) 2 . In 
the present work a new approach has been developed permitting more precise calculations of the energy. In particular, 
the terms of the order of (D / J) would give us a contribution of the Dzyaloshinskii-Moriya interaction into the energy 
of the cubic anisotropy responsible for (i) the ordering of the spiral axes along selected crystallographic directions, 
e.g. (Ill) in MnSi, in the absence of the external magnetic field; (ii) the orientation of the triangle Skyrmion lattice 
in the A-phase observed in these crystals 
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Appendix A 

In this work we often need to take sums over twelve bonds ij of components of the vectors Tj and Tj . It is convenient 
to calculate the sums with a computer, but in simple cases we can also use the formulae 



V Ny N M M y = N N y N A ( v-4 M M y M z 

Z-^{ij} ix iy iz jx jy jz \A^i— 1 ix iy iz J I Z-/j=l jx jy jz 

_y-4 N x +M x N y 
1 ix iy 



(Al) 



4 



V" n x n v Nz _ ) 4, N x (mod2) = N y (mod 2) = N z (mod2), 
0, m other cases. 



1=1 



The vectors Tj, Tj, (n - Tj), (r, + Tj), [r, i x Tj], i>n,ij, ^n.ij ( n — 1,2,3) change with the use 
of symmetry transformations of the point group 23, when the index ij passes trough 12 possible values 
{12, 13, 14, 21, 23, 23, 31, 32, 34, 41, 42, 43}. Therefore we can use the formulae 

^>y=0, (A3) 
{«} 



a ^»b l3 [3 = 4(ai2 • b 12 )S al 3, (A4) 

{ij} 

where ay and by are the vectors from the above-listed set. Thereby the following summations can be easily made: 



2 1 = 12, (A5) 

{ij} 



J2(Ti±Tj) = 0, 

{ij} 

X Tj] =0, 

{ij} 



(A6a) 
(A6b) 
(A6c) 
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/ ^ \ , i ' j }o>\' i \ > j } p 


(A7a) 






? ITi zb T^,-vItV X Twlfl = 

/ ,\ 't — *J/Q^L * J IP ' 








7 (TV — T« 1^ IT* — T» ) = 32(5^ 








^(Tj + Tj) a (n + Tj)p = 16d a p, 


(A7d) 


{y} 




^[Ti X Tj] Q [Ti X = 32<5 Q/3 , 


(A7e) 






y^CTjetTj/g + TjaTifj) — ~85 a p. 


(A7f) 



Using Eqs. (IA7a V (|A7b ) we can easily calculate sums containing products of the vectors h nyi j, H n ,ij (n = 1,2,3): 



(A8a) 

V"(t» - Tj) a b 2 ,ijp = 32 ( | - x J (5 Q/3 , (A8b) 
V^Ti - Tj) a h,ijp = 32 ( -i - x J <5 Q/3 , (A8c) 



/ — Tj) a D nii j@ = 8D n+ S al 3, (A9) 



HP = 4£>, 2 1 (Sq/3. (A10) 



If the parameter x is related with w by Eq. ([5]), then 



^ bi tija bi t ijff = ( - + 2w 2 )(5 Q(3 , (Alia) 
{«} 

2J b 2 ,i ja b2,ij0 =(--4w + 2w 2 )S af3 , (Allb) 

b3.ijab3.ijp = h + Aw + ^w 2 )5 a p, (Allc) 



and 



^ D nt i ja b n j jf3 = (-2wD n+ + 4D„ 42 ■ b^ 12 )£ a /3, ( Al2 ) 

{ij} 



or 
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D i,VctKw = (-2w£>i+ - D x - + 2D ly )8 a p, (A13a) 

to} 

D 2t i ja b2,ijp = (-2^2+ + 2D 2+ + £> 2 - + 2D 2y )8 a p, (A13b) 

to} 

D 3,ijab 3 ,m = {~2wD 3+ - 2D 3+ + + 2D 3y )5 a p, (A13c) 

to} 



where 12 is the vector b n ,i2 in the case of w = 0. 
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Figures 
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FIG. 1: (color online), (a) The ID twisted magnetic structure; the spins rotate in the plane perpendicular to the propagation 
vector, (b) The canting in the absence of a twist; the lattice is divided into two sublattices (cyan and sand-coloured arrows) 
with the spins tilted by a constant angle, (c) The canting and the twist; the canting of the sublattices leads to the combination 
of two helices with the same propagation vector and a constant phase shift. 
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FIG. 2: (color online). Nearest magnetic neighbourhood of a manganese atom in MnSi. The central atom (magenta) belongs to 
the first manganese sublattice, its neighbors are in 2nd (orange), 3rd (green) and 4th (blue) ones. The shortest bonds are given 
by solid lines, the bonds with the atoms from the 2nd and 3rd magnetic shells are shown by dashes and dots, correspondingly, 
(a) The 3-fold axis [111] is vertical. The horizontal lines designate equidistant atomic planes containing the manganese atoms 
of 2nd, 3rd and 4th sublattices. (b) The projection on the plane perpendicular to the axis [111]. Two atomic triangles from 
(a), the upper and the lowest, belonging to different shells are situated just one above the other. 
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FIG. 3: (color online), (a) The ambiguity of the transition from discrete to continuous model: two different smooth functions 
having the equal values in a discrete set of points, (b) The gain from the arbitrary choice of the parameter x. The initial curve 
(orange) is plotted with the use of the function values given in the discrete set of points with the coordinates x and x within 
the cells of a ID crystal. When the points are shifted to the center of the cells without a change of function values, the view 
of the curve changes (blue). 




FIG. 4: The areas of the ferromagnetic (dark) and the antiferromagnetic (light) order in the space of the exchange interaction 
parameters j± — (J2 ± Js)/ Ji- The case of J\ > is presented. When Ji < 0, the areas exchange their colors. 
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(b)Pfi X*! 13 1/2 = 




FIG. 5: (color online). Periodical ID chain of magnetic atoms with two spins within a unit cell. The spins can rotate in a plane, 
and their directions are fully determined by the only variable cp. The spin interaction is described by the classic Hamiltonian 
(|55|) . (a) The bonds of the atom 1 of the nth unit cell and the parameters defining the interactions between the spins. A 
chirality is given by the difference (C2 — Ci). (b) When J2 = 0, the chirality is hidden; the lattice is divided onto two sublattices 
with the ferromagnetic ordering of spins and a canting angle /3. (c) When a ferromagnetic interaction with the atoms of the 
alternate kind from the second neighboring cells is included (J2 > 0), the chirality becomes apparent in the magnetic helix with 
the angle S between neighboring unit cells; besides, the angle a between the atoms 1 and 2 whithin a cell remains the same 
over the chain. 



